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We investigate the 5=1/2 quantum spin antiferromagnets on the triangular and Kagome lattices 
in magnetic field, using the numerical exact diagonalization. Particularly we focus on an anomalous 
magnetization behavior of each system at 1/3 of the saturation magnetization. The critical exponent 
analyses suggest that it is a conventional magnetization plateau on the triangular lattice, while an 
unconventional phenomenon, called the magnetization ramp, on the Kagome lattice. 
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The S = 1/2 triangular and Kagomei lattice anti- 
ferromagnets have attracted a lot of interests as typ- 
ical frustrated systems. Most theoretical studies in- 
dicated that the former system has the three sublat- 
tice long-range order— ~—, while the latter is disordered 
in the ground stated—. Experimental studies to ob- 
serve a novel spin liquid phase have been accelarated 
since discoveries of several realistic materials; the organic 
compound k-(BEDT-TTF) 2 Cu 2 (CN) 3 for the triangu- 
lar lattice^, the herbertsmithit o 18 ' 19 , the volborthit o 20 ' 21 
and the vesignieite 2 ^ for the Kagome lattice. Since the 
quantum Monte Carlo simulation and the DMRG cal- 
culation are useless for these systems, the numerical ex- 
act diagonalization is one of the best numerical method 
for them. The numerical diagonalization studies sug- 
gested that both systems have the 1/3 magnetization 
plateau 2 ^— , although the classical spin systems have no 
plateau on both lattices in the ground stat o 28 i 29 . (The 
thermal or quantum fluctuations induce a plateau in the 
semiclassical case, because 1/3 is just a critical point be- 
tween two dirrerent spin structures.) In our recent nu- 
merical diagonalization study on the 5 = 1/2 Kagome 
lattice antiferromagnet up to N = 36, the calculated 
field derivatives reveal an anomalous behavior at 1/3 of 
the saturation magnetization 30 . Namely, the field deriva- 
tive is diversing at the low-field side of the critical field 
H c , while almost zero at the high-field side. This critical 
behavior is quite different from conventional magnetiza- 
tion plateaux in two-dimensional systems where the field 
derivative is finite at both sides of H c . To distinguish 
such an anomalous property at the 1/3 magnetization of 
the Kagome lattice from conventional plateaux, we called 
it a "magnetization ramp" . However, its mechanism is 
still an open problem. In this paper, to clarify such an 
unconventional behavior around the 1/3 magnetization of 
the S =1/2 Kagome lattice antiferromagnet, comparing 
with the triangular one, we applied the numerical diag- 
onalization for both systems up to N = 39 which is the 
largest cluster at present. In addition we estimated the 



critical exponent 5 by the finite-size scaling proposed by 
the previous worksa, to investigate the quantum critical 
behevior more quantitatively. 

Now we examine the magnetization process of the S — 
1/2 triangular and kagome lattice antiferromagnets. The 
Hamiltonian is given by 



ri — T~Lq + H2 



(1) 
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n a =^s l -s„ n z = -nY J S], (2) 

where means all the nearest neighbor pairs on each 
lattice. Throughout we use the unit such that gfiB = 1- 
For TV-site systems, the lowest energy of Hq in the sub- 
space where J^j &j — M (the macroscopic magnetization 
is rn = M/M s , where M s denotes the saturation of the 
magnetization, namely M s = NS for the spin-S* system. 
) is denoted as E(N, M). We restrict us to the rhombic 
cluster under the periodic boundary condition to keep 
the 120° rotational symmetry for a systematic finite-size 
scaling. Using the numerical exact diagonalization, we 
have calculated all the values of E(N, M) available for the 
rhombic clusters with N =9, 12, 21, 27, 36 and 39, to ob- 
tain the ground state magnetization curves. (The largest 
dimension of the N = 39 system is 68,923,264,410. To 
treat such huge matrices in computers, we have carried 
out parallel calculations using the MPI-parallelized code 
which is originally developed in the previous work^l.) 
The magnetization curves of the (a) triangular and (b) 
Kagome lattice antiferromagnets are shown in Fig. [1] 
for N —27, 36 and 39. They indicate plateau-like be- 
haviors at m — 1/3 of both systems, but the Kagome 
lattice exhibits an anomalous feature; the step length in- 
creases with decreasing H towards m = 1/3, different 
from conventional magnetization plateau like the trian- 
gular lattice. In order to clarify a difference between the 
triangular and Kagome lattice, we calculated the field 
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FIG. 1: Magnetization curves of the (a) triangular and (b) 
Kagome lattice antifferomagnets for N =27, 36, 39. 



derivative \ defined the form 

_! _ E(N, M + 1) - 2E(N, M) + E(N, M + 1) 
X ~ 1/M. • ( j 

The derivative x °f the (a) triangular and (b) Kagome 
lattice systems are shown in Fig. [2] for N =27, 36 and 
39. The derivative x °f the triangular system is finite at 
both edges of the 1/3 plateau-like behavior, as well as 
conventional two-dimensional systems. In contrast, the 
Kagome system exhibits a quite different feature between 
the lower and higher field sides of m = 1/3; x is diverg- 
ing at the lower side like a plateau in one-dimensional 
systems, while is very small (possibly zero) at the higher 
one. The present calculation for N = 39 more strongly 
supports a ramp-like behavior predicted by our previous 
work. 

The critical exponent 6 defined by the form 



\H - H, 



1 1/5 



(4) 



is an important index to specify the universality class of 
the field induced quantum phase transition. The previ- 
ous theoretical works indicated 6 = 2 for some typical 
one-dimensional gapped system s 33 i 34 , while 6 = 1 for 
two-dimensional systems 35 . In order to investigate the 
quantum critical behavior at m = 1/3 of the triangu- 
lar and Kagome lattice antiferromagnets more quantita- 
tively, we estimate <5 by the finite-size scaling developed 
by the previous work 32 . Although it was proposed for 



0.6 



0.4 



I 1 

- (a) 


1 1 1 1 1 


1 


- 






P' ' 

- Ap'' 




(J [ 
-~o~ 


f 




— -A--- 


1 1 


s 27 

1 1 1 1 1 


~-a~- 



0.5 

M/M s 



1.5 



0.5 



1 1 1 1 1 1 1 

: (b) ? 




if 




aq 




*** . 

i i i i l i i 





0.5 1 

M/M s 

FIG. 2: Field derivatives of the (a) triangular and (b) Kagome 
lattice antifferomagnets for TV =27, 36, 39. 



one-dimensional systems, it can be easily generalized for 
two dimensions. We assume the asymptotic form of the 
size dependence of the energy as 



±-E(N,M)~e(m) + C(m)j^ 



(N^oo), (5) 



where e(m) is the bulk energy and the second term de- 
scribes the leading size correction. We also assume that 
C(m) is an analytic function of m. The lowest and high- 
est magnetic field corresponding to m = 1/3 in the ther- 
modynamic limit are defined as H c - and H c+1 respec- 
tively, as the form 

N N 
E{N,-±l)-E(N,-)^±H c± (N-too). (6) 

In order to consider the critical magnetization behaviors 
for m < 1/3 and m > 1/3 independently, we define the 
critical exponents 6 C - and 6 C + by the forms 



|m-i|~|#-tfc±| 1/5± . 



(7) 



If we define the quantities /+ (N) and /_ (N) by the forms 
f±(N) 

ee ±[E(N, j ± 2) + E(N, j) - 2E{N, j ± 1)],(8) 
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the asymptotic forms of them are expected to be 
1 „, 1 



f±(N) 



N 5 ± 



+ o(^tt) (n 



(9) 



as far as we assume the form ([5]). Thus the exponents 
6- and 8+ can be estimated from the slope of the In f±- 
\nN plot, respectively, under the condition 9 > 6± — 1. 
In order to avoid an oscilation of the finite-size correc- 
tion due to the cluster shape dependence, we just use the 
rhombic clusters under the periodic boundary condition 
with N =9, 12, 21, 27, 36, and 39. The plots of ln/± 
versus In N for the (a) triangular and (b) Kagome lattice 
antiferromagnets are shown in Fig. [31 Figure [3J (a) sug- 
gests that the calculated points are well fitted to a line 
for each of /_ and /+ in the case of the triangular lat- 
tice. Thus applying the standard least square fitting to 
lines (Dashed and long-dashed lines are used to obtain 8 + 
and 5-, respectively.) for all the available system sizes; 
N =9, 12, 21, 27, 36, 39 (N = 9 cannot be used for S-), 
5- and 8+ are estimated as follows: 



<5_ = 1.00 ±0.17, 



0.89 ±0.15, 



for the triangular lattice. The errors are estimated from 
the deviation of points from the fitted lines. It would 
be reasonable to conclude <5_ = 8+ = 1 at m = 1/3 
of the triangular lattice antiferromagnet, as expected for 
conventional magnetization plateau in two dimensions. 
On the other hand, Fig. [3] (b) indicates quite different 
feature of the Kagome lattice antiferromagnet. The same 
least square fitting yields the following estimates: 



6_ = 1.92 ±0.99, 



0.56 ±0.15, 



for the Kagome lattice antiferromagnet. Exponent 8- 
has a large error because the line fitting is not good. It 
does not converge with respect to the system size well, 
but seems to still increasing with N. The same line fitting 
to the points for N =27, 36 and 39 yields the estimation 
<5_ = 4.59 ± 0.25. Thus we can just conclude <5 > 2 
at most. It means that the diversing behavior of the 
field derivative at 7J C _ is stronger than one-dimensional 
systems. It is led to two possibilities. One is a jump 
(a first-order transition) in the magnetization curve. A 
magnetization jump which also appears near the satura- 
tion was proved^. The other is an anomalous continu- 
ous transition. A similar phenomenon was reported in 
the metal-insulator transition of the Hubbard chain with 
next-nearest neighbor hopping^. In comparison with <5_ , 
(5+ is more conclusive, because the fitting error is much 
smaller. According to the above result of the line fitting, 
we conclude 8+ = 0.6 ± 0.2. Thus the field derivative x 
should be zero at the higher field side of H c+ , because 8+ 
is smaller than unity. It also justifies a property of the 
magnetization ramp. 

Finally, we consider whether a flat part of the magne- 
tization curve at m = 1/3 exists or not for the triangular 
and Kagome lattice antiferromagnets. Namely, we exam- 
ine whether each system has no plateau (H c - = H c+ ) or 
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FIG. 3: ln(/) is plotted versus ln(JV) for (a) the triangular and 
(b) Kagome lattice antiferromagnets, respecitively. Solid cir- 
cles and squares are useful to estimate the critical exponents 
8+ and 8—, respectively. 



a finite plateau (H c - ^ H c+ ) at m = 1/3 in the thermo- 
dynamic limit. We evaluate the length of the flat part 
H c + — H c - corresponding to the plateau width of the 
finite-size clusters with N =9, 12, 21, 27, 36 and 39 for 
both systems. If the system has a gapless excitation like 
a spin wave from some ordered states, the low-lying en- 
ergy spectrum is expected to be proportional to the wave 
vector k in the long wave lenght limit. Thus the excita- 
tion energy gap of the finite-size systems should have 
the asymptotic form ~ 1/iV 1 / 2 in two-dimensional gap- 
less systems. On the other hand, in gapped systems the 
gap is expected to converge to the thermodynamic limit 
with exponentially decaying (faster than 1/7V 1 / 2 ) finite- 
size correction, as the system size increases. Thus if the 
extrapolation by fitting the gap versus 1/N 1 / 2 leads to 
a finite gap in the thermodynamic limit, it would be a 
strong evidence to confirm the gapped ground state. The 
length of a flat part H c+ — H c - is plotted versus l/N 1 / 2 
in Fig. 21 where open triangles and solid circles are for 
the triangular and Kagome lattice antiferromagnets, re- 
spectively. The least square fitting to a line leads to 
the following results: H c+ - H c _ = 0.47 ± 0.28 for the 
triangular lattice and H c+ — H c - = —0.32 ± 0.35 for 
the Kagome lattice. Obviously we can conclude that the 
triangular lattice antiferromagnet has the 1/3 magneti- 
zation plateau. In contrast, the result for the Kagome 
lattice suggests that it possibly has a single critical field 
H c = H c - = H c+ . However, it is difficult to exclude a 
finite magnetization plateau, because of a large error of 
the extrapolation. Note that any other plateaux are diffi- 
cult to investigate by the present method, because fewer 
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FIG. 4: Plateau width H c+ - H c - is plotted versus 1/A 1/2 . 
Open triangles and solid squares are for the triangular and 
Kagome lattice antiferromagnets, respectively. Fitted lines 
are used for the extrapolation to the thermodynamic limit. 

system sizes can be available for m ^ 1/3. 

In the recent magnetization measurement^ on a can- 
didate of the Kagome lattice antiferromagnet Volborthite 
several step-like behaviors were observed, but it has not 
reached m = 1/3 yet. The same measurement is sitll 
going on to observe an anomaly at m = 1/3, which is ex- 
pected to be about 60T. It would be interesitng to detect 
some unconventional features. 

In summary, we have investigated critical magnetiza- 



tion behaviors at m = 1/3 for the S = 1/2 triangu- 
lar and Kagome lattice quantum antiferromagnets, us- 
ing the numerical exact diagonalization of rhombic clus- 
ters up to JV = 39. The triangular lattice is revealed 
to have the critical exponets 5— = 5+ = 1 and a finite 
plateau, which are consistent with a conventional mag- 
netization plateau in two-dimensional systems. On the 
other hand, the Kagome lattice is revealed to exhibit un- 
conventional critical properties; 5_ < 1 < <5 + , namely 
the field derivative \ is diverging at the lower field side, 
while zero at the higher one of a possibly single critical 
field H c = H c - = H c+ . The conclusion supports the 
magnetization ramp behavior at m — 1/3 of the Kagome 
lattice antiferromagnet. 
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